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Theoretical  and  Experimental 
Studies  of  HF  Ducted  Propagation 


1.  INTRODUCTION 

The  propagation  of  high  frequency  (HF)  radio  energy  in  so-called  "non- 
classical"  modes,  that  is,  in  modes  other  than  the  conventional  earth-ionosphere 
hops,  has  recently  attracted  a  great  deal  of  attention.  This  report  deals  with  two 
such  propagation  modes,  the  chordal  and  ducted,  often  grouped  together  under  the 
generic  term  "ducting."  Figure  1  illustrates  the  two  different  modes  schemat¬ 
ically;  note  that  the  chordal  mode  is  a  special  case  of  a  ducted  mode.  Numerical 
calculations  and  computer  propagation  simulations  are  described.  New  data  on 
the  details  of  the  vertical  profile  of  the  lower  ionosphere  are  incorporated  into 
some  of  these  calculations.  The  question  of  how  HF  energy  can  be  coupled  into 
elevated  ducts  from  a  ground  based  source  is  examined  in  detail.  Absorption 
losses  in  ducts  are  estimated  using  specially  adapted  raytracing  procedures. 

A  number  of  different  but  related  experiments  are  described  in  which  attempts 
have  been  made  to  measure  some  of  the  properties  of  ducted  propagation  modes. 

In  some  of  these,  natural  ionospheric  gradients  were  used  to  inject  energy  into 
ducts,  while  in  others  a  h:gh  power  ionospheric  heater /modifier  was  used  as  a 
means  of  mode  coupling. 

In  1976,  the  concept  emerged  in  the  Propagation  Branch  of  the  Electromag¬ 
netic  Sciences  Division  to  examine  the  possibility  of  using  low -loss  ionospheric 

(Received  for  publication  21  November  1980) 
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Figure  1.  Schematic  Display  of  Chordal  and  Ducted  Mode 


ducting  channels  for  propagating  electromagnetic  signals  to  large  distances.  The 
prediction  for  the  availability  of  such  naturally  occurring  ducts,  based  on  idealized 
global  ionospheric  conditions  had  already  been  established. 1  Viewing  the  iono¬ 
sphere  as  a  refractive  medium,  one  observes  that  favorable  horizontal  ionization 

gradients  are  needed  for  rays  to  enter  an  elevated  duct  from  a  ground-based 
2 

source.  In  practice,  such  favorable  gradients  are  available  at  certain  hours  of 
the  day  depending  on  location  and  their  estimated  local  occurrence  rate  is  about 
12  percent.  Favorable  gradients  in  the  general  vicinity  of  a  transmitter  may 
lead  to  the  capture  of  a  sufficient  portion  of  the  electromagnetic  signal  in  a 
naturally  occurring  duct.  To  receive  such  a  signal  on  the  ground  at  great  dis¬ 
tances  from  the  source  would  require  either  favorable  positive  gradients  associ¬ 
ated  with  horizontal  ionization  enhancement  in  the  direction  of  propagation  that 


1.  Tushentsova,  I.  A.,  Fishchuk,  D.  I. ,  Tzedilina,  Ye.  Ye.  (1975)  Investigation 

of  the  global  properties  of  ionospheric  wave  ducts.  2.  Geomaen.  Aeron. 
15(l):62-66.  - - 

2.  Toman,  K. ,  Miller,  D.  C.  (1977)  Computational  study  of  long-range  high- 

frequency  ionospheric  ducting.  Radio  Science  12(3):467-476. 
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could  bend  a  ray  enough  to  reach  the  ground,  or  ionospheric  irregularities  that 
might  scatter  a  substantial  portion  of  the  electromagnetic  signal  to  a  receiver 
at  the  ground. 

Prevailing  ionospheric  ducting  channels  may  not  always  be  accessible  to 

transmissions  originating  from  the  ground.  Local  modification  of  the  ionosphere 

by  radio  frequency  heating  creates,  however,  ionospheric  irregularities  that 

could,  by  a  scattering  process,  deflect  radio  energy  into  an  existing,  elevated 
3 

ionospheric  duct.  Conversely,  electromagnetic  signals  already  trapped  in  such 
a  long-range  duct  could  be  ejected  from  the  duct  toward  the  ground  by  natural 
and  artificially  induced  ionospheric  irregularities. 

While  in  certain  applications  ionospheric  modification  may  eventually  be 
required  near  both  ends  of  a  propagation  path,  the  study  described  here  was 
conceived  to  examine  the  feasibility  of  successfully  using  an  available  radio 
frequency  heater  to  eject,  with  artificially  created  ionospheric  irregularities, 
ducted  signals  to  the  ground.  For  the  path  from  Adelaide,  Australia  to  Los 
Alamos,  advantage  was  taken  of  favorable  ionospheric  gradients  near  the  trans¬ 
mitter  site  at  sunset.  Conditions  were  enhanced  by  favorable  gradients  associ¬ 
ated  with  the  equatorial  anomaly  in  the  F-region  ionosphere. 


2.  GENERAL  THEORETICAL  CONSIDERATIONS 

In  order  to  estimate  the  characteristics,  availability,  and  behavior  of  iono¬ 
spheric  ducting  channels,  several  techniques  are  available,  each  having  its  own 
advantages  and  limitations.  These  techniques  can  be  distinguished  as  1)  Numer¬ 
ical  raytracing  through  model  ionospheres  considered  representative  of  condi¬ 
tions  occurring  in  nature;  2)  qr-diagram  obtained  for  a  representative  ionospheric 
profile  of  a  spherically  symmetric  medium  without  magnetic  field  and  without 
collisions;  3)  potential  field  analogy  for  a  representative  ionospheric  profile  of 
a  spherically  symmetric  medium  and  its  generalization  to  apply  to  an  inhomo¬ 
geneous  medium.  Departures  from  spherical  symmetry  generally  pose  no  prob¬ 
lems  to  numerical  raytracing  but  tend  to  invalidate  the  use  of  the  potential  field 
approach  when  inhomogeneities  become  severe.  The  qr-diagram  remains  strictly 
applicable  to  a  spherically  symmetric  medium  only. 

It  is  believed  that  the  limitation  imposed  on  the  validity  of  raytracing  is  less 
stringent  than  the  constraint  of  the  medium's  spherical  symmetry.  Over  the 
short  distances  within  which  sudden  changes  in  the  refractive  index  would  restrict 


3.  Gurevich,  A.V. ,  Tzedilina,  Ye.  Ye.  (1976)  Trapping  of  radiation  in  the 
ionospheric  duct  during  scattering  on  artificial  inhomogeneities, 
Geomagn.  Aeron.  15(6):713-715. 
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the  use  of  raytracing,  specification  of  the  medium  is  not  easily  accomplished  in 
practice.  When  a  propagation  medium  is  specified  for  the  purpose  of  raytracing, 
care  is  taken  in  applying  mathematical  algorithms  and  profile  specification  func¬ 
tions  for  multiple  layers  to  assure  that  the  refractive  index  varies  slowly  in  all 
directions. 

2.1  Potential  Reid  Analogy 

The  propagation  of  radio  waves  in  ionospheric  ducts  may  be  treated  in  a 

4 

manner  analogous  to  the  trapping  of  nuclear  particles  in  a  potential  field  well. 

Let  U(h)  =  -e(h)  describe  the  field  intensity  as  a  function  of  altitude,  h, 
where 


e(h)  =  1 


->) 


2h 

ir 


and 

u  =  wave  frequency 

uQ  =  ionospheric  plasma  frequency  ~  ^Ng(h) 

Ng(h)  =  vertical  profile  of  electron  density 

R  =  earth  radius 

o 

If  o(h)  is  the  local  elevation  angle  of  the  ray  direction  at  height  h,  then  by 
Snell's  law: 

2 

e(h)  cos  a(h)  =  constant 

For  energy  trapped  in  an  E-F  duct,  with  reflection  at  h^,  then  o(h1)  =  0,  and  the 
maximum  value  of  a  in  the  duct  is  given  by 


4.  Gurevich,  A.  V.  (1971)  Effect  of  nonlinearity  on  the  generation  of  circum¬ 
terrestrial  signals,  Geomagn.  Aeron.  11(6):810-817. 
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An  adiabatic  invariant  can  be  similarly  defined,  for  a  frequency  f,  as 


««  =  - /  Ie(h,  f)  -  e  ,  (h,  f)] 1/2  dh 


The  value  of  the  adiabatic  invariant  is  proportional  to  the  energy  density  in  the 
propagating  mode  at  the  frequency  in  question. 

Figure  2  illustrates  the  form  of  the  function  e(h)  for  a  typical  vertical  profile. 
Ng(h).  In  this  example,  energy  may  be  trapped  in  two  altitude  regions  corre¬ 
sponding  to  two  different  ranges  of  the  potential  function.  These  two  ducting 
modes  are  illustrated  in  Figure  1,  in  which  radio  energy  is  shown  respectively 
trapped  between  the  E  and  F  layers  of  the  ionosphere  and  propagating  via  a 
glancing  or  chordal  mode  on  the  underside  of  the  F-layer. 

Figure  3  shows  the  average  ionospheric  electron  density  profile  in  mid¬ 
latitudes  as  measured  by  a  direct  ion  probe  on  board  the  satellite  S3-1.  This 
profile,  as  indicated,  is  typical  of  the  latitude  range  45°  to  60°  at  low  magnetic 
activity  level  (Kp  =  0-2)  and  for  solar  zenith  angles  (70°  <  x  <  90°).  Note  the 


ELECTRON 
DENSITY-N#(h ) 


\/Vi 

A  R- 


CHORDAL  MODE 


E-F  DUCT 


Figure  2.  Schematic  Display  of  Potential  Function  for  Two-layer  Ionosphere 


R.  Philbrick,  Private  Communication, 
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Figure  3.  Mid-latitude  Electron  Density  Height  Profile  Averaged 
from  Direct  Ion  Probe  Measurements  on  Board  the  Satellite  S3-1 


well-defined  "valley"  at  about  180  km,  between  the  E  and  F  regions,  which  can 
be  expected  to  support  a  ducted  mode.  Figure  4  shows  the  value  of  e(h)  plotted 
as  a  function  of  altitude  for  the  profile  of  Figure  3,  for  several  values  of  radio 
frequency  between  10  and  30  MHz.  The  potential  well  at  170  km  altitude  can  be 
seen  to  decrease  in  intensity,  as  expected,  as  the  frequency  increases. 

While  the  description  of  radio  wave  ducting  in  terms  of  a  potential  well  is  a 
useful  one,  there  is  this  limitation:  that  the  potential  function  depends  on  one 
spatial  variable  only,  and  thus  only  a  spherically  stratified  ionosphere  can  be 
considered.  This  simplification  excludes  the  treatment  of  coupling  of  energy 
between  a  ducted  mode  and  a  ground-based  transmitter  or  receiver  (in  Figure  1, 
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Figure  4.  Potential  Functions  for  the  Ionosphere  of  Figure  3 
That  Display  the  Changes  With  Operating  Frequency 


the  source  of  energy  is  shown  located  in  the  duct).  An  advantage  of  the  potential 
well  analogy  is  its  simplicity,  which  permits  the  convenient  estimation,  for  exam¬ 
ple,  of  the  global  ducting  properties  of  the  ionosphere  from  maps  of  predicted 
ionospheric  critical  frequencies. 

2.2  Illustration  of  Potential  Field  Technique* 

In  order  to  illustrate  the  value  of  the  potential  field  technique  in  estimating 
the  likelihood  of  ionospheric  ducting  over  large  distances,  an  example  is  pre¬ 
sented  here.  A  transmitter  is  postulated  at  a  mid-latitude  location  (28°N;  183°E) 
with  a  beam  oriented  45°  West  of  North  (azimuth  315°).  The  ionosphere  was 
modelled  along  this  path  using  the  parameters  of  the  CCIR  ionospheric  model. 
Models  were  constructed  for  an  equinox  month  (March)  at  00,  06,  12,  18  UT  for 
Sunspot  Numbers  (SSN)  of  50  (medium)  and  150  (high).  Figure  5  shows,  for  exam¬ 
ple,  the  UT  dependence  of  the  maximum  and  minimum  critical  frequencies  (f  F„) 

O  £ 

along  the  propagation  path.  Figure  6  shows  the  variation  with  range  along  the 

direction  of  propagation  of  the  height  profile  of  plasma  density,  at  UT  =  00, 

SSN  =  150.  Range  is  expressed  in  earth  centered  degrees  (180°  ~  22,  000  km) 

-3  -5 

and  plasma  density  In  (cm  X  10  ).  At  UT  =  0,  the  transmitter  is  close  to  local 

noon;  sunset  occurs  at  a  range  of  about  70°,  and  the  equatorial  anomaly  is  seen 
roughly  centered  at  about  130°  in  the  nighttime  ionosphere.  Figures  7  and  8 
show  contours  of  fQFg  as  functions  of  range  and  UT  for  Sunspot  Numbers  150  and 
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Figure  5.  Computed  Time  Variation  of 
the  Maximum  and  Minimum  of  the  F- 
region  Critical  Frequencies  ((^2)  Along 
a  Particular  Propagation  Path  for  Medium 
and  High  Sunspot  Number 


50,  respectively.  Also  indicated  on  these  plots  are  the  regions  of  day  and  night, 

the  latter  shown  shaded.  The  altitude  variation  of  plasma  frequency  along  the 

direction  of  propagation  is  shown  in  an  isometric  format  in  Figure  9  for  00  UT  at 

SSN  =  150.  The  maximum  values  of  the  three  layers  (f  E,  f  F,  and  f  F„)  are 

J  o  ’  o  1  o  2 

indicated  by  heavy  solid  lines. 

The  potential  function  corresponding  to  Figure  9  is  shown  in  Figures  10,  11, 
and  12  for  frequencies  20,  22,  and  24  MHz,  respectively.  From  the  transmitter 
out  to  a  range  of  7,  200  km,  in  the  daytime  ionosphere,  classical  hop  type  E  and 
F  modes  are  in  evidence  as  well  as  E-F  ducted  energy,  depending  on  the  elevation 
angle  and  frequency.  At  night  (for  ranges  greater  than  7,  200  km)  only  classical 


16 


CONTOURS  OF  N#X  K>S(CM"3) 


MARCH 
SSN  ISO 
UTOO 


90 

RANGE  (DEG) 


Figure  6.  Computed  Electron  Density  Height  Contours  Along  a  Particular 
Propagation  Path 


F-mode  and  ducted  F  (that  is,  chordal)  propagation  is  possible.  In  Figure  10, 
corresponding  to  a  frequency  of  20  MRz,  the  boundaries  of  the  various  modes  are 
indicated  by  heavy  solid  lines;  progressing  from  the  lowest  altitudes  upwards, 
these  lines  indicate  (1)  lower  boundary  of  chordal  E  mode,  (2)  upper  boundary 
of  chordal  E,  or  lower  boundary  of  E-F  ducted  mode,  (3)  upper  boundary  of 
E-F  ducting,  and  (4)  upper  boundary  of  chordal  F.  For  example,  in  Figure  10, 
an  E-F  duct  exists  out  to  range  7,  200  km,  while  a  chordal  E  mode  is  possible 
between  3,  300  km  and  6,  600  km  ranges.  At  night,  only  ground  hop  F  modes 
exist  except  between  8,  300  km  and  11,  600  km,  when  a  chordal  F  mode  may  prop¬ 
agate.  It  should  be  noted,  in  this  example  at  least,  that  the  E-F  duct  potential 
well  is  very  shallow  and  in  fact  is  scarcely  discernible  to  the  eye.  It  is  also  clear 
that  a  chordal  E  mode,  while  theoretically  possible  in  the  daytime,  would  suffer 
very  high  absorption  losses  in  practice. 
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Figure  7.  Computed  Contours  of  F-region  Critical 
Frequencies  (fQF2)  Spanned  by  Distance  and  Time  (UT) 
for  a  Particular  Propagation  Path  (SSN  =  150) 

The  altitude  range  over  which  ducted  propagation  is  possible  is  determined 
by  the  maximum  valve  (that  is,  the  smallest  negative  value)  of  the  potential 
function  for  the  frequency  of  interest.  Figure  13  shows  a  plot  of  this  maximum 
value  for  the  E-layer  and  F-layer  boundaries,  for  selected  frequencies  (16,  18 
and  20  MHz)  for  UT  =  00,  SSN  =  150.  In  the  free  space  region  below  the  iono¬ 
sphere,  chordal  modes  are  only  possible  when  t  <  -1.  It  is  apparent  that  decreas 
ing  the  frequency  reduces  the  range  over  which  chordal  modes  may  propagate. 

In  fact  at  night,  in  the  example  Illustrated,  no  chordal  F  propagation  may  exist 
for  frequencies  less  than  about  16.  5  MHz. 


SSN50 


UT  (HR) 

Figure  8.  Computed  Contours  of  F-region  Critical 
Frequencies  (fQFo)  Spanned  by  Distance  and  Time  (UT) 
for  a  Particular  Propagation  Path  (SSN  =  50) 

An  estimate  of  the  energy  density,  which  the  various  ducted  modes  are  ca¬ 
pable  of  supporting  along  the  path,  is  obtained  by  computing  the  adiabatic  invariant 
defined  above.  Figure  14  shows  1(f)  plotted  along  the  direction  of  propagation  for 
a  frequency  of  20  MHz.  Note  the  different  ordinate  scales  for  the  F  (chordal)  and 
E/E-F  modes,  respectively,  illustrating  the  fact  that  the  F-Mode  can  carry  about 
20  times  more  energy  than  the  chordal-E  or  E-F  ducted  modes. 


PLASMA  FREQUENCY 


Figure  9.  Model  Plasma  Frequency  for  Three -layer  Ionosphere  Varying 
With  Height  Above  and  Range  Along  a  Particular  Propagation  Path 
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Figure  11.  Isometric  Display  of  Potential  Function  for  Ionosphere  of 
Figure  9  Using  Operating  Frequency  of  22  MHz 
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Figure  12.  Isometric  Display  of  Potential  Function  for  Ionosphere  of 
Figure  9  Using  Operating  Frequency  of  24  MHz 


Figure  14.  Plots  of  the  Adiabatic  Invariant  for  E,  EF  and  F 
Mode  Varying  Along  the  Propagation  Path.  (The  ordinate 
scale  for  the  F  mode  differs  from  that  of  the  E,  EF  modes. ) 


2.3  Coupling  to  Duct  Mode  Propagation 

In  the  case  of  a  smooth,  spherically  stratified  ionosphere,  it  is  geometrically 
impossible  for  radio  energy  to  be  coupled  from  a  ground  based  source  to  a  ducted 
propagation  mode;  chords  of  the  spherical  ionospheric  cavity  which  satisfy  the 
reflection  condition  and  intersect  the  earth's  surface  at  one  point  (the  source) 
will  always  intersect  the  earth  at  other  points  as  well.  In  order  to  couple  energy 
from  a  ground  based  source  into  a  ducted  mode,  some  perturbation  of  the  ray 
trajectory  is  required  at  the  coupling  point.  This  perturbation  may  take  one  of 
several  forms:  (1)  refraction  due  to  a  suitable  ionospheric  gradient;  (2)  scatter¬ 
ing  from  natural  small-scale  irregularities  (for  example,  auroral  or  equatorial 
irregularities);  (3)  scattering  from  artificially  produced  ionospheric  irregular¬ 
ities  (for  example,  those  induced  by  RF  heating,  or  those  accompanying  release  of 
suitable  chemicals  or  plasma  in  the  ionosphere,  or  those  produced  by  electro¬ 
static  acceleration  of  ionospheric  electrons).  Other  possibilities  exist,  such  as 
a  combination  of  ionospheric  gradients  and  scatter  from  meteor  trails  in  the 
lower  ionosphere.  ® 

5.  Vever,  A.S. ,  Danilova,  T.  P. ,  Shlionsky,  A.G.  (1978)  Possibility  of  the 

trapping  of  radio  waves,  scattered  by  meteor  trails,  by  ionospheric  ducts, 
Geomagn.  Aeron.  18(3):306-308. 


25 


2.4  Raytracing  Through  Natural  Gradients 

In  general,  a  decreasing  electron  density  gradient  in  the  direction  of  propa¬ 
gation  is  required  to  couple  energy  from  a  ground  based  source  into  a  ducted 

mode  (conversely,  an  increasing  gradient  is  required  for  coupling  from  a  ducted 

2 

mode  to  a  ground  based  receiver).  Suitable  gradients  exist  under  commonly 
occurring  situations,  such  as  crossing  the  solar  terminator  (day-night  boundary 
in  the  ionosphere),  propagation  through  the  polar  and  equatorial  regions,  and 
traveling  ionospheric  disturbances.  The  frequency  range  and  the  range  of  eleva¬ 
tion  angles,  over  which  trapping  in  ducts  may  take  place,  are  functions  of  the 
magnitude  of  the  gradient  and  of  the  location  of  the  irregularity  relative  to  the 
source  on  the  ground. 

Computations  have  been  performed,  using  models  of  the  ionosphere,  together 
with  numerical  raytracing  procedures,  to  establish  the  range  of  frequencies  and 
elevation  angles  over  which  ducting  can  occur  in  typical  cases.  An  example  of 
such  a  computation  is  presented  here  in  which  the  influence  of  strong  gradients 
near  sunset  is  considered. 

A  transmitter  was  assumed  to  be  directing  radio  energy  at  low  elevation 
angles,  from  a  location  at  latitude  32°  in  a  northeasterly  direction.  The  iono¬ 
sphere  was  modeled  along  the  direction  of  propagation  at  1415,  1515,  1615  and 
1715  LT  at  the  transmitter,  for  the  month  of  October,  at  sunspot  number  100. 
Figure  15  shows  the  location  of  the  evening  solar  terminator,  at  these  four 
times,  relative  to  the  transmitter  and  direction  of  propagation.  It  is  seen  that  at 
1415  the  entire  ionospheric  propagation  path  is  in  sunlight,  while  at  1715  it  is  in 
darkness.  At  1515  and  1615  l.Tf  large  gradients  are  modeled  along  the  propaga¬ 
tion  path. 


Figure  15.  Progression  of  Solar  Terminator 
Toward  Transmitter 


Three  dimensional  numerical  raytracing  was  performed  at  a  frequency  of 
30  MHz  in  the  four  modeled  ionospheres,  with  the  results  depicted  in  Figure  16 
where  the  top  panel  corresponds  to  1415  LT  and  time  progresses  downwards  to 
1715  in  the  bottom  panel.  Each  panel  shows  an  elevation  cut  and  an  azimuthal 
cut  through  the  fan  of  rays  from  0°  to  10°  in  steps  of  1°  in  elevation.  The  lower 
part  of  each  panel  shows  the  elevation  cross  section  relative  to  the  surface  of  the 
earth,  while  the  upper  part  shows  the  azimuthal  cross  section  relative  to  the 
propagation  direction.  Range  is  indicated  by  the  scale  graduated  in  units  of 
500  km  (the  first  500  km  is  subdivided  into  100  km  steps). 


Figure  16.  Sequence  of  Three-dimensional  Numerical  Raytracings 
Through  Model  Ionosphere  at  Four  Local  Times  for  a  Fixed  Trans¬ 
mitter 


The  daytime  and  nighttime  ionospheres  display  essentially  classical  hop  type 
propagation  (note  that  the  ray  plots  are  terminated  when  the  ray  intersect  the 
earth).  At  1715  (night),  however,  one  ray  is  apparently  ducted  and  a  second  just 
grazes  the  earth.  The  intermediate  ionospheres  support  injection  of  the  entire 
transmitted  rays  into  chordal  and  ducted  modes,  as  a  result  of  the  large  negative 
electron  density  gradients  along  the  direction  of  propagation. 


2.5  Scattering  from  Irregularitiee 


The  phenomenon  of  scattering  of  radio  waves  from  ionospheric  small-scale 
irregularities  is  essentially  the  same  for  both  naturally  occurring  and  artificially 
produced  irregularities.  The  only  practical  differences  lie  in  the  possibly  differ¬ 
ent  scale-size  spectra  and  different  spatial  extents  involved  in  the  two  casea 
Figure  17  illustrates  the  geometry  involved  when  radio  waves  of  power  density 
Pq  are  incident  on  elongated  irregularities  at  angle  x-  The  power  scattered  into 
an  element  of  solid  angle,  dfj,  at  angle  to  the  irregularity  axis,  from  an  incre¬ 
mental  element  of  scatterer  volume,  dV,  is 

dPs  =  Po  0<O)  dn  dV 

where  6  is  the  scattering  angle,  that  is,  the  angle  between  incident  and  scattered 
vectors  and  <r(0)  is  the  "volumetric  scatter  cross  section,"  that  is,  the  equivalent 
point  target  cross  section  per  unit  volume  of  scatterers.  In  the  ionosphere,  the 
scatterers  are  elongated  along  the  direction  of  the  local  magnetic  field;  the 
volumetric  cross  section  is  a  strong  function  of  the  scatter  angle  0. 

The  equivalent  point  target  cross  section  of  a  finite  volume  Vs  of  scatterers 
is 


as(0)  -  / 

V 


a(0)  dV 


0 

Erukhimov  et  al  have  given  the  formula  for  the  volumetric,  cross  section  of 
long,  thin  plasma  irregularities  having  a  Gaussian  electron  density  distribution, 
for  the  general  case  of  arbitrary  angles  of  incidence  and  scatter.  Taking  their 
results  and  making  minor  changes  in  symbols,  we  obtain 


o(0,  y)  = 


kV 

16\^r 


<(£e)2)  •  exp  -  (krjf  sin  |-)2 


,  2  ,  1 
(cos  y  +  ~j 

V 


sin2  y) 


where 

k  =  2n/\ 

A  =  wavelength 


6.  Erukhimov,  L,  M. ,  Matyugin,  S.  N.,  Uryadov,  V.  P.  (1975)  Radio  wave 
propagation  in  ionospheric  wave  channels.  Radtophysics  &  Quantum 
Electronics  18(9):958-963.  ' 
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i  =  scale  size  transverse  to  the  magnetic  field 
2 

((Ae)  )  =  mean  square  fluctuation  of  dielectric  constant  corresponding  to 
scale  size  I 


r)  =  ratio  of  irregularity  dimensions  parallel  and  transverse  to  the 
magnetic  field 


cos  7  =  CP-S-1  --c°£  X 
2  sin  ^ 


Note  that 


<->2>  ■  $  w 

where 

uQ  =  local  plasma  frequency 
u  =  2irc/\  =  frequency  in  radians 
c  =  velocity  of  light 

N  =  electron  density  of  background  plasma 
AN  =  rms  electron  density  fluctuation 
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The  expression  for  volumetric  cross  section  may  be  written  in  the  useful  form 


cr(0,  4>,  X) 


4  .3 

k 

16\fjr 


T 

w 


exp  <-  k“T 


1) 


(cos  0  -  COS  x) 


2 


+  sin 


There  are  two  important  special  cases  of  the  formula  for  volumetric  cross  sec¬ 
tion.  (1)  Scattering  cone:  When  0  =  x,  the  first  term  in  the  exponential  factor 
is  zero,  corresponding  to  a  maximum  in  cross  section.  This  is  the  so-called 
specular  reflection  condition,  and  is  satisfied  on  the  surface  of  a  cone  whose  axis 
is  the  axis  of  the  irregularity.  (2)  Forward  scatter:  When  0  =  0,  the  second 
factor  In  the  exponential  term  is  zero,  corresponding  to  a  maximum  in  cross 
section. 
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Figure  17.  Scattering  Geometry  for  Reradiation  of  Incident  Electromagnetic 
Waves  by  Elongated  Irregularities 


An  interesting  special  case  is  that  in  which  the  incident  and  scattered  vectors 
are  both  in  the  magnetic  meridian  plane,  in  which  case 

0  =  0  -  X 

and  the  volumetric  cross  section  is  a  function  of  ^  and  x  alone. 

Figures  18  and  19  illustrate  the  variation  in  reflected  power  (normalized  to 
unit  incident  power)  as  a  function  of  scatter  angle  0  for  ranges  of  values  of  the 
transverse  irregularity  dimension  and  anisotropy.  Figure  20  illustrates  the 
aspect  sensitivity  of  the  scatter  for  the  same  ranges  of  parameters,  as  a  function 
of  wavelength,  close  to  orthogonality.  [  Note  that  the  curve  for  (/  =  100,  r?  =  20) 
coincides  almost  exactly  with  that  for  (f  =  300,  n  =  10).]  It  is  apparent  that  the 
aspect  sensitivity  of  the  scatter  increases  as  irregularity  size  and  anisotropy 
increase  (for  a  given  wavelength)  or  wavelength  decreases  (for  given  irregularity 
dimensions). 
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Figure  20.  Computed  Wavelength  Dependence  of  Aspect 
Sensitivity  Near  Orthogonality  (0  -»  0) 


2.  5.  1  SCATTERING  AT  HIGH  FREQUENCIES 

At  frequencies  well  above  wo>  the  plasma  frequency,  the  refraction  of  incident 
and  scattered  waves  is  small,  and  straight  line  propagation  geometry  is  a  reason¬ 
ably  accurate  approximation.  If  attention  is  also  confined,  for  simplicity,  to  . 
propagation  in  the  magnetic  meridian  plane,  <//  is  the  local  magnetic  inclination 
for  ducted  waves  traveling  horizontally. 
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2.  5.2  SCATTER  GEOMETRY  WITH  RADIO  WAVE  REFRACTION 

The  field-aligned  scatter,  which  can  act  as  a  means  of  coupling  energy  in  and 
out  of  ionospheric  ducts,  has  thus  far  been  considered  only  in  the  high  frequency 
limit,  that  is,  with  consideration  of  radio  refraction  neglected.  In  fact,  for  radio 
frequencies  in  the  HF  band,  considerable  refraction  is  encountered  in  the  radio 
path  through  the  ionosphere  below  the  scattering  region,  especially  at  the  low 
end  of  the  band  and  in  the  daytime.  Thus  the  straight-line  propagation  geometry 
must  be  modified  to  reflect  the  actual  propagation  conditions,  which  is  a  task  of 
considerable  complexity. 

Computations  have  been  performed  in  which  the  ionosphere  was  modeled  in 
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three  dimensions,  using  a  general  purpose  model,  and  rays  were  traced  appro¬ 
priately  through  this  model  over  a  range  of  frequencies.  A  ducted  mode  was 
assumed  to  be  propagating  towards  the  scattering  point  from  one  of  several 
azimuths  (0°  to  360°  in  45°  increments);  after  specular  scattering  and  refracted 
propagation  through  the  lower  ionosphere,  the  points  of  intersection  of  the  rays 
with  the  earth's  surface  were  computed.  In  order  to  permit  some  deviation  from 
a  rigid  requirement  for  horizontal  propagation  in  the  ducted  wave,  the  angle  of 
incidence  upon  the  scattering  point  was  assumed  to  take  the  values  0°,  ±5°  with 
respect  to  the  horizontal  at  that  location.  This  permitted  some  estimate  to  be 
made  of  the  ground  locus  of  intersection  within  a  10  degree  range  of  angles  for  the 
ducted  wave  vector.  (See  Section  3.  2,  Figure  33.) 

2.6  Propagation  Loaa  Evaluation 

In  order  to  estimate  the  absorption  losses  for  radio  waves  trapped  in  iono¬ 
spheric  ducts,  a  propagation  analysis  was  conducted  in  which  a  numerical  iono¬ 
spheric  model  was  constructed  and  rays  subsequently  traced  in  this  model.  The 

modeling  and  raytracing  procedure  which  was  used  permitted  the  computation  of 

7 

absorption  due  to  both  electron-neutral  and  electron-ion  collisions.  Radio  waves 
propagating  in  the  F-region  are  influenced  relatively  little  by  electron-neutral 
collisions,  due  to  the  relatively  low  concentration  of  neutral  atoms  at  F-region 
altitudes;  the  major  contribution  to  absorption  at  these  altitudes  comes  from 
electron-ion  collisions. 

A  source  of  radio  energy,  in  the  frequency  range  10  to  21  MHz  was  assumed 
to  exist  at  several  altitudes  in  the  ionosphere,  and  propagation  to  a  range  of 
10,  000  km  was  considered.  The  elevation  angle  of  the  initial  ray  vector  at  the 


7.  Rush,  C.M.,  Elkins,  T.J.  (1975)  An  assessment  of  the  magnitude  of  the 
F-region  absorption  on  HF  radio  waves  using  realistic  electron  density 
and  collision  frequency  models,  ITU  Telecommunication  Journal.  Geneva. 
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source  was  varied  over  a  wide  range  and  this  angle,  together  with  the  ionospheric 
structure  along  the  path,  determined  the  propagation  mode.  Figure  21  shows  the 
results  of  the  attenuation  calculations  for  four  values  of  the  source  altitude  be¬ 
tween  110  and  170  km.  Also  illustrated  are  two  possible  propagation  modes  —  one 
classical  hop-type  and  one  ducted  (chordal).  The  general  trend  in  the  dependence 
of  attenuation  on  frequency  reflects  the  approximate  inverse-square  dependence 
which  is  well-known  from  classical  theory.  There  is  also  a  pronounced  inverse 
dependence  upon  the  altitude  of  the  source,  which  reflects  the  larger  proportions 
of  the  propagation  paths  in  the  lower  regions  of  the  ionosphere  for  lower  source 
altitudes.  For  the  170-km  source  altitude,  a  group  of  points  is  found  to  exhibit 
extremely  low  absorption  in  the  frequency  range  15  to  20  MHz,  as  indicated  by 
the  dashed  contour.  These  correspond  to  ducted  modes,  for  which  the  absorption 
losses  are  in  the  range  2  ~  4  dB  for  the  10,  000-km  path. 


Figure  21.  Sample  Illustration  of  Propagation  Loss,  Computed  for  Different 
Transmitter  Heights,  Over  a  10,  000-km  Propagation  Path 


3.  GROUND-BASED  PROPAGATION  EXPERIMENTS 
3.1  Round-the-Wortd  (RTW)  Propagation 

High  frequency  waves  which  travel  completely  round  the  world  (so-called 
RTW  propagation)  may  be  taken  as  evidence  of  the  existence  of  ionospheric  duct¬ 
ing  over  appreciable  portions  of  the  path;  otherwise  lower  ionospheric  absorption 
in  the  sunlit  hemisphere  would  have  reduced  their  signal  strength  below  the  level 
of  detectability.  A  related  indication  of  ducting  is  the  propagation  of  significantly 
large  energy  along  the  long  portion  of  the  great  circle  joining  two  relatively  close 

points  (so-called  "back-path"  propagation).  Round-the-world  propagation  has 
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been  studied  for  many  years,  but  a  consistent  explanation  of  the  phenomenon  has 

g 

only* recently  emerged.  Based  on  theoretical  and  experimental  considerations, 

it  is  apparent  that  the  optimum  conditions  for  RTW  propagation  are  along  a  great 
circle ‘which  is  close  to  the  solar  terminator  (the  twilight  line).  Generally  speak¬ 
ing,  the  greater  the  angular  separation  between  the  propagation  path  and  the 
optimum  RTW  great  circle,  which  itself  is  displaced  15°  to  20°  from  the  termi¬ 
nator,  the  less  favorable  are  conditions  for  RTW  propagation. 

A  number  of  RTW  and  related  experiments  have  been  conducted  to  investigate 
the  ducting  phenomenon.  An  FM/CW  transmitter,  operating  at  1-kW  power  in  the 
frequency  range  6  to  30  MHz  was  coupled  to  a  horizontally  polarized  log  periodic 
antenna  (gain  — 14  dB  relative  to  isotropic)  directed  northwards  from  Ava, 

New  York  (lat.  45°).  The  receiving  antenna,  located  at  nearby  Verona,  New  York 
consisted  of  an  array  of  Beverage  elements  having  a  mid-band  gain  of  about  12  dB, 
directed  south.  Swept  frequency  ionograms  were  recorded  around  the  clock  at 
20  min  intervals,  with  RTW  propagation  being  frequently  observed.  Figure  22 
shows  an  example  of  such  an  ionogram;  note  the  diffuse  nature  of  the  trace  and 
the  bifurcation  at  the  high  frequency  end,  a  characteristic  feature  of  these  RTW 
measurements.  The  maximum  and  minimum  observed  frequencies  were  extracted 
from  a  sample  eight  hour  period  of  ionograms  and  are  plotted  in  Figure  23  to¬ 
gether  with  the  observed  group  delay.  Also  shown  are  the  time  of  sunrise  at  the 
approximate  latitude  of  the  ionospheric  entry  point  (55°)  and  at  the  equator  of  the 
longitude  of  the  transmitter/ receiver  pair.  Note  the  approximate  symmetry  of 

the  period  of  RTW  observation  about  12  l)T,  that  is,  along  the  optimum  RTW 

g 

great  circle.  The  group  path,  measured  at  the  time  of  arrival  of  the  earliest 
signal  returns  (the  leading  edge)  shows  a  maximum  variation  of  about  0.  5  msec 


8.  Toman.  K.  (1979)  High-frequency  ionospheric  ducting  -  A  review.  Radio 

Science  14(3):447-453. 

9.  Golyan,  S.  F.  (1975)  Optimal  conditions  for  extra-long-range  short-wave 

propagation,  Radiophysics  A  Quantum  Electronics  18(9):1014-1022. 
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Figure  22.  Round-the-world  Echo  of  FM/CW  Signal  Radiated  Northward  From 
Ava,  NY  and  Received  at  Verona,  NY  About  137.5  msec  later 


(150  km)  over  the  observing  period.  Since  the  longer  group  paths  are  associated 
with  the  highest  values  of  the  minimum  observed  frequency,  it  is  presumed  that 
the  propagation  path  involved  more  ground  reflections,  resulting  in  higher 
absorption,  during  the  period  ~  11  to  14  UT  than  in  the  earlier  part  of  the  ob¬ 
served  RTW  period. 

A  second  type  of  measurement  was  made  in  which  transmissions  from  Ava 
were  directed  both  towards  the  east  and  towards  the  west,  with  reception  at 
Los  Alamos,  New  Mexico.  The  Platteville  heating  facility,  operated  by  the 
Institute  for  Telecommunications  Sciences,  Boulder,  Colorado,  was  used  to 
modify  the  ionosphere  with  its  ten  200  kW  transmitters,  feeding  a  circular  array 
of  broadband  dipoles  that  can  tune  from  5  to  10  MHz.  With  an  antenna  gain  of 
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Figure  23.  Observed  Maximum  and  Minimum  Frequencies 
in  Megahertz  of  Round-the-world  Echoes  and  Their  Group 
Delays  in  Milliseconds 


20  dBi  and  a  half-power  beamwidth  of  16  degrees  at  7.  5  MHz  the  heater  could 
modify  a  400  X  400  X  20  (altitude)  km  volume.  It  was  of  interest  to  determine 
whether  the  Platteville  heater  was  operating  during  this  period,  and  the  intent 
was  to  determine  whether  or  not  radio  energy  could  be  scattered  from  an  iono¬ 
spheric  duct  forming  a  part  of  the  RTW  propagation  path. 

Figures  24  and  25  illustrate  two  examples  of  quasi-RTW  propagation.  In  the 
first  example,  Ava  is  transmitting  towards  the  east  so  that  the  energy  is  propa¬ 
gating  along  the  "back-path";  in  the  second  example,  Ava  is  transmitting  west, 
and  the  received  radio  energy  has  propagated  slightly  longer  than  one  transit 
around  the  earth.  These  two  examples  were  selected  to  show  two  different  types 
of  RTW— the  first  exhibiting  substantial  multipath  (that  is,  several  relatively 
discrete  modes)  whereas  the  second  is  more  intense  and  has  a  sharply  defined 
leading  edge  with  diffuse  trailing  edge.  During  both  of  the  measurements  illus¬ 
trated  here,  the  ionospheric  heater /modifier  was  off. 
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Figure  24.  Round-the-world  Echo  of  FM/CW  Signal  Radiated  Eastward 
from  Ava,  NY  and  Received  at  Los  Alamos,  New  Mexico 


Figures  26  to  30  show  examples  of  data  sets  in  which  the  heater  was  cycled 
on  and  off  at  5-min  intervals  or  longer  (indicated  by  the  binary  vertical  markers 
along  the  abscissas).  Again  plotted  are  the  upper  and  lower  observed  RTW  fre¬ 
quencies  and  the  group  delay.  Most  of  these  data  were  collected  near  the 
autumnal  equinox  and  Figure  31  shows  the  distribution  of  observed  RTW  periods, 
grouping  around  sunset  (~00  UT)  and  sunrise  (~  14  UT).  Figure  30  shows 
summer  data,  with  RTW  echoes  centered  around  sunset  ( —  02  UT).  Also  shown 
in  Figure  30,  for  comparison  purposes,  is  the  maximum  observed  frequency 
(MUF)  propagating  over  the  direct  path  from  Ava  to  Los  Alamos,  indicating  that 
the  maximum  RTW  frequency  was  only  1  to  2  MHz  higher  than  the  classical  MUF 
for  most  of  the  observing  period.  It  appears  from  the  data  in  Figures  26  to  30 
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Figure  25.  Round -the -world  Echo  of  FM/CW  Signal  Radiated  Westward 
From  Ava,  NY  and  Received  at  Los  Alamos,  New  Mexico 
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that  the  presence  of  heater-generated  ionospneric  irregularities  has  no  discern¬ 
ible  effect  on  the  range  of  frequencies  propagating  in  the  RTW  signal.  Qualitative 
examination  of  the  data  reveals  that  the  heater  also  has  no  apparent  effect  on  the 
intensity  of  the  RTW  signals. 

The  observation  that  heating  of  the  ionosphere  had  no  apparent  effect  on  the 
RTW  signal  has  at  least  two  possible  interpretations:  (1)  The  scatter  cross 
section  of  the  heater -produced  irregularity  cloud  was  sufficiently  small  that  any 
additional  energy  scattered  from  the  RTW  signal  was  below  the  noise  level;  this 
is  not  unreasonable  since  the  measured  RTW  levels  were  very  close  to  noise 
level.  (2)  The  irregularity  cloud  was  too  high  in  the  ionosphere  with  the  result 
that  the  heater-induced  scatterers  were  outside  the  range  of  altitudes  in  which 
the  ducted  signal  was  propagating. 
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Figure  28.  Sep  20,  1978,  Ava  -  Los  Alamos,  Observed  Maximum  and 
Minimum  Frequencies  in  Megahertz  of  Round-the-world  Echoes  and 
Their  Group  Delays  in  Milliseconds 


Figure  29.  Sep  22,  1978,  Ava  -  Los  Alamos,  Observed  Maximum  and 
Minimum  Frequencies  in  Megahertz  of  Round-the-world  Echoes  and 
Their  Group  Delays  in  Milliseconds 
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figure  30.  June  5,  1979,  Ava  -  Los  Alamos,  Observed  Maximum  and 
Minimum  Frequencies  in  Megahertz  of  Round-the-world  Echoes  and 
?heir  Group  Delays  in  Milliseconds 
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Figure  31.  Preferred  Times  of  Occurrence  of  Round-the-world  Echoes 
Between  Ava,  New  York  and  Los  Alamos,  New  Mexico  from  18  to  22 
Sep  1978 
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The  second  possibility  is  regarded  as  being  the  more  likely  one  since  the 
lower  frequency  limit  of  the  heater  (5  MHz)  was  close  to  the  maximum  F-region 
plasma  frequency  for  much  of  the  nighttime  observing  periods.  Thus  heating 
was  taking  place  mostly  above  about  250  km,  and  RF  energy  propagating  in  the 
chordal  mode  (the  dominant  ducting  mode  at  night)  can  be  expected  to  travel  at 
a  lower  altitude  than  this  value. 

3.2  Long-range  Propagation 

Several  ducted -propagation  experiments  have  been  conducted  between  1978 
and  1980  which  involved  the  ionospheric  radio-frequency  heating  facility  at 
Platteville,  Colorado,  operated  by  the  Institute  for  Telecommunications  Sciences, 
Boulder,  Colorado.  This  radio -frequency  heating  facility  was  used  to  create  in 
the  F-region  a  heater  cloud  whose  fine  structure  contained  field-aligned  ioniza¬ 
tion  irregularities.  These  irregularities  were  expected  to  scatter  incoming  high- 
frequency  signals,  originating  from  remote  transmitting  sites,  to  receiving  sites 
whose  general  location  relative  to  Platteville  was  deduced  from  the  geometry 
provided  by  ionization  irregularities. 

Characteristics  of  the  geometry  were:  (1)  Irregularities  are  constrained  to 
be  aligned  with  the  orientation  of  the  local  geomagnetic  field,  (2)  irregularities 
occur  at  heights  that  depend  on  the  height  distribution  of  the  ambient  ionosphere 
plasma  frequency  and  the  heater  frequency,  (3)  Direction  of  the  incident  radiation 
originating  from  remote  ground-based  transmitters. 

Transmitting  and  receiving  sites  used  in  the  experiments,  along  with  surface 
distances  and  propagation  delay  between  stations  are  illustrated  in  Figure  32. 
Table  1  gives  the  location  of  participating  stations,  transmitter  power,  antenna 
gain,  and  range  of  operating  frequencies. 

The  sounding  transmissions  from  the  Stanford  Research  Institute,  Inter¬ 
national  facility  at  Lost  Hills,  California,  from  the  RADC  facility  at  Ava,  New 
York,  and  from  the  Australian  site  at  Adelaide  employed  a  linear  swept -frequency 
continuous  wave  (SFCW)  system.  These  transmissions  were  directed  toward 
Platteville,  Colorado,  where  the  RF  heater  transmitter  was  operated  by  the 
Institute  for  Telecommunications  Sciences.  The  most  notable  effect  of  the  heat¬ 
ing  process  is  the  creation  of  ionospheric  irregularities  in  the  F  region  by  the 
heater  beam.  The  artificially  created  irregularities  within  the  heated  volume  are 
aligned  in  the  direction  of  the  earth's  magnetic  field.  A  receiving  station  located 
at  Los  Banos,  California  obtained,  for  the  Lost  Hills  transmission,  backscatter 
information  on  the  cross  section  of  the  irregularity  volume  over  Platteville,  as 
well  as  ground  backscatter  by  way  of  ionospheric  reflections.  The  receiving 
stations  at  Los  Alamos,  New  Mexico  and  Alamosa,  Colorado,  were  selected  on 
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ADELAIDE 

Figure  32.  Layout  for  High-Frequency  Propagation  Experiment  Utilizing  the 
Ionospheric  Heater.  (Units  for  surface  distance  and  transmission  delay  are 
kilometer  and  millisecond;  solid  points  represent  transmitters;  open  points 
represent  receiving  stations. ) 


Table  1.  Participating  Transmitting  and  Receiving  Stations 


Location 

Latitude 

Longitude 

Transmitter 

Power 

(kW) 

Antenna 

Gain 

(dBi) 

Frequency 

(MHz) 

Platteville, 

Colorado 

40.  2°N 

104. 7°W 

1500 

20 

5-10 

Adelaide, 

Australia 

34. 85°N 

138. 66°E 

10 

20 

6-30 

Los  Banos, 
California 

37.  06°N 

120.  83°W 

- 

28 

6-30 

Lost  Hills, 
California 

35.  8°N 

120°W 

20 

18 

6-28 

Alamosa, 

Colorado 

37. 43°N 

105.  86°W 

■IBBB 

10 

6-30 

Los  Alamos, 
New  Mexico 

35,  76°N 

106,  23°W 

■  ■ 

10 

6-30 

Los  Lunas, 
New  Mexico 

34.82° 

106.  74°W 

16 

6-30 

ANew  York 

43°N 

75.  5°W 

10 

16 

6-30 

the  basis  of  computations  of  the  family  of  scattering  cones,  of  field-aligned  irreg¬ 
ularities,  that  intersect  the  surface  of  the  earth  at  loci  which  depend  on  the  chosen 
geomagnetic  field  model,  arrival  angle  of  the  electromagnetic  wave  relative  to 
magnetic  azimuth  and  elevation  angle  incident  at  the  irregularity,  the  height  of 
the  irregularity,  and  sounding  frequency  because  of  frequency-dependent  refrac¬ 
tion  effects  in  the  ionosphere. 

Figure  33  illustrates  the  results  of  sample  calculations  of  loci  resulting  from 
the  intersection  of  scattering  cones,  associated  with  field-aligned  irregularities, 
with  the  surface  of  the  earth,  refraction  in  the  ionosphere  included.  The  condi¬ 
tions  derived  from  an  ionospheric  model  are  typical  of  autumn,  noon,  sunspot 
number  150.  The  origin  of  the  scattering  cones  are  placed  at  a  height  of  200  km. 
At  that  height  above  Platteville  (P),  Colorado  (lat.  40.  2°N,  long.  104.  7°W;  see 
Table  1)  electromagnetic  signals  were  allowed  to  arrive  at  elevation  angles  of 
0°,  +5°,  -5°.  In  Figure  33,  this  angle  was  +5°  above  the  horizontal  at  a  height 
of  200  km  over  Platteville. 

Signal  frequencies  of  6,  12,  18,  24  and  30  MHz  were  used  to  compute  rays 
that  are  specularly  reflected  from  field-aligned  irregularities,  each  frequency 
resulting  in  a  family  of  discrete  rays  per  angle  of  incidence  that  defines  a 
"refracted"  scattering  cone.  The  azimuth  of  rays  were  0,  45,  90,  135,  180° 
relative  to  the  orientation  of  the  magnetic  meridian  (see  Figure  34)  with  east  or 
west  azimuth  providing  identical  loci  for  intersection  between  scattering  cone 
and  surface  of  the  earth.  While  the  magnetic  field  was  used  to  define  the  scat¬ 
tering  geometry,  it  was  neglected  in  the  computations  for  the  rays.  Rays  were 
computed  for  all  frequencies  and  azimuths  cited  earlier.  Not  all  loci  were 
plotted,  however,  in  Figure  33.  For  certain  combinations  of  frequency  and 
azimuth,  it  was  found  that  the  scattering  cones  did  not  intersect  the  surface  of 
the  earth,  for  example,  24-  and  30-MHz  signals  at  azimuth  135°M.  For  the 
same  frequencies  arriving  at  an  azimuth  of  45°M,  however,  the  signals  are 
scattered  to  the  ground  south  of  Los  Alamos. 

For  ionospheric  conditions  typical  of  autumn,  local  midnight,  two  intersec¬ 
tions  of  the  scatter  cone  including  refraction  are  illustrated  for  18  MHz  and 
azimuths  90°M,  135°M.  For  the  conditions  chosen  in  Figure  33  it  is  seen  that 
the  locus  of  intersection  for  18  MHz  (135°M)  differs  significantly  from  that  shown 
in  Figure  35  for  day  versus  night,  200  versus  250  km  altitude,  and  assumed 
incident  elevation  angle  +5°  versus  0°.  For  local  noon,  autumn  and  a  200-km 
altitude  at  which  signals  at  12,  18,  24  and  30  MHz  coming  from  an  azimuth  of 
135°M  are  assumed  to  be  incident  on  field -aligned  irregularities.  Figure  36 
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Figure  33.  Loci  of  Scattering  Cones,  Including  Refraction,  Intersect¬ 
ing  the  Surface  of  the  Earth.  Scatter  altitude:  200  km.  Irregularities 
are  field  aligned,  located  above  Plattevllle,  Incident  ray  arrives  at 
irregularity  height  with  elevation  angle  of  5°  from  above  the  horizon¬ 
tal.  Azimuth  (AZ)  of  incidence  is  counted  relative  to  magnetic  meridian. 
This  azimuth  angle  is  zero  if  the  ray  arrives  from  the  north  in  the  plane 
of  the  magnetic  meridian 
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Figure  34.  The  Azimuthal  Configuration  Over  Platteville, 
Colorado,  Pertaining  to  the  Scattering  Loci,  the  Magnetic 
Meridian  and  the  Direction  Toward  Adelaide,  Australia 
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Figure  35.  Loci  of  Scattering  Cones  with  Refraction  for  f  =  18  MHz 
and  2  Azimuths  Illustrating  the  Location  of  the  Loci  North  and  South 
of  Los  Alamos 
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Figure  36.  Effect  of  Elevation  Angle  of  Ray  Arriving  at  Irregularity 
Height  with  Elevation  Angle  of  +5°  (from  above)  and  -5°  (from  below) 
With  Respect  to  the  Horizontal 


illustrates  the  effect  of  elevation  angle  at  the  points  of  incidence  above  Platteville. 
The  solid  contours  refer  to  arrival  angle  above  the  horizontal  of  +5°,  the  dotted 
contours  to  -5°.  For  signals  assumed  to  arrive  from  the  south-west  (135°M),  a 
negative  elevation  angle  of  -5°  displaces  the  intersection  loci  farther  south  from 
Platteville  compared  with  a  positive  angle  of  +5°.  Higher  frequencies  tend  to 
come  farther  south  because  they  encounter  reduced  refraction. 

Figures  33,  34  and  36  were  based  on  simulation  data  derived  from  an  iono¬ 
spheric  model  that  provided  the  geometry  of  the  scattering  contours.  The 
results  illustrate  some  characteristics  of  contours  obtained  by  the  choice  of 
model  ionospheres  representing  mean  conditions.  The  choice  of  initial  conditions 
for  electromagnetic  signals  arriving  at  the  field-aligned  columns  of  the  scatter¬ 
ing  volume  resulted  in  an  infinitesimal  scattering  cone  that  provided  the  initial 
conditions  for  computing  the  rays  of  which  some  reached  the  surface  of  the  earth, 
while  others  did  not. 

The  scattering  geometry  which  is  in  principle  characterized  by  the  magnetic 
field  of  the  earth  is  practically  time  invariant.  In  practice,  ionospheric  condi¬ 
tions  will  generally  be  different  from  the  mean  model  used  for  computations. 

These  changing  conditions  will  influence  the  choice  of  heater  frequency  and  thus 
the  altitude  at  which  field  aligned  irregularities  are  produced.  Its  consequence 
will  be  a  set  of  different  intersection  loci  originating  now  from  a  particular 
scattering  cone  located  at  some  altitude  with  the  landing  of  rays  being  subject  to 
the  prevailing  refractive  influences  of  the  ionosphere. 

3.  2.  1  MEASUREMENT  CAMPAIGN  OF  1978 

During  July  and  September  1978  measurements  were  made  of  the  HF  induced 
backscatter  over  Colorado**  with  an  over-the-horizon  radar  at  Lost  Hills,  and  the 
Wide  Aperture  Research  Facility  (WARF)  at  Los  Banos,  both  in  California.  The 
radar  operated  as  a  linear  swept-frequency  continuous  wave  (SFCW)  backscatter 
sounder  at  frequencies  from  6  to  30  MHz.  The  ionospheric  modification  was  pro¬ 
duced  by  the  Platteville  heater  facility  using  5  to  10  MHz  as  modifier  frequencies. 
The  most  notable  ionospheric  effect  of  radio  frequency  heating  is  the  creation  of 
irregularities  in  the  heater  beam  which  are  aligned  in  the  direction  of  the  earth's 
magnetic  field.  Narrow -band  FM/CW  transmissions  (about  25-kHz  bandwidth) 
were  used  to  measure  "cloud"  echo  strength  resulting  from  the  radar  illumination 
of  the  irregularity  volume  and  observe  cloud  cross-section  changes  caused  by  step¬ 
wise  variations  in  the  heater  power  levels,  A  difference  in  transmitted  power  of 
about  3.  5  dB,  resulting  from  cycling  the  RF  heating  power  between  1.  20  and 
0.  54  MW  in  2  min  intervals,  resulted  in  a  variation  of  echo  power  by  about  2  dB 

11.  Showen,  R.L.,  Zavoll,  W.B.  (1979)  Platteville  Heating  Results:  HF  Radar 
Observations.  SRI  International,  Tech.  “Report  No.  ft. 
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indicating  a  saturation  with  increasing  HF  power  of  the  process  causing  the 
scatter. 

At  the  WARF  receiving  site  in  Los  Banos,  the  use  of  a  hybrid  spectrum 
analyzer  provided  real-time  wide-sweep  backscatter  soundings  to  aid  in 
optimizing  the  selection  of  transmissions  and  recording  parameter^.  The 
receivers  were  operated  with  automatic  gain  control  to  accommodate  the  very 
large  variation  in  received  signal  strength  typical  of  wide-sweep  soundings. 

To  obtain  azimuthal  profiles  of  the  "cloud's"  cross  section,  three  of  the  WARF's 
half-degree  receive  beams  were  repetitively  stepped  across  a  series  of  contin¬ 
uous  bearings  as  the  sounding  frequency  was  varied.  A  WARF  backscatter  sound¬ 
ing  ionogram  showing  "direct"  scatter  from  the  heated  region  is  illustrated  in 
Figure  37.  Calibration  markers  were  superimposed  at  intervals  of  1  msec  in  time 
delay  and  1  MHz  in  frequency  corresponding  to  10  seconds  for  the  sweep  rate  of 
100  kHz  per  second.  Identifiable  signatures  in  the  sounding  include  the  F-layer 
land -backscatter  returns  and  direct  echoes  from  the  heated  volume  (cloud)  over 
Platteville.  The  2 -MHz  periodic  structure  seen  in  the  "cloud"  echo  is  the  result 
of  stepping  the  receive  antenna  beam  in  azimuth  every  2.  5  seconds  corresponding 
to  a  2-MHz  interval.  The  echoes  from  the  "cloud"  seem  to  extend  in  frequency 
over  the  16  to  30  MHz  band.  The  mechanism  of  standard  field-aligned  scatter 
responsible  for  the  cloud  echoes  led  to  absolute  cross-section  estimates  of  10® 
square  meters.  Although  more  information  about  the  scattering  cross-section 
is  available,  the  cross-section  estimates  generally  agree  with  earlier  results. 12 

The  computation  of  the  apparent  cross-section  a  which  includes  an  excess 
loss  term  that  cannot  be  controlled  or  directly  measured  by  the  radar,  was  con¬ 
veniently  introduced.  The  equation  for  this  quantity  is 


=  r  =  S 


4*R2 

PtGt 


4  it 
G  X2 


4*TR 


where 

<xc  =  radar  backscatter  cross  section  of  the  "cloud" 

L  =  excess  loss 

S  =  received  signal  strength 

R  =  slant  range 

Pt  =  transmitted  power 


12.  Flaler,  P.A.  (1974)  Field-aligned  scattering  from  a  heated  region  of  the 
ionosphere  -  Observations  at  HF  and  VHF,  Radio  Science  9(ll):923-940. 
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=  receive  antenna  gain 
X  =  wavelength  of  sounding  signal 


a 
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apparent  cross-section 
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Figure  37.  Wide  Aperture  Research  Facility  (WARF)  Backscatter  Sounding  (Shows 
Direct  Scatter  From  the  Heated  Volume  Above  Platteville  and  the  Effects  of 
Azimuthal  Scanning  on  the  Delay  Spread) 


3.  2.  2  MEASUREMENT  CAMPAIGN  OF  1979 
3.  2.  2.  1  Backscatter  Ionograms 

The  second  set  of  coordinated  propagation  experiments  was  conducted  in  April 
and  June  1979  with  RADC,  Stanford  Research  Institute,  International  (SRI), 
Institute  for  Telecommunication  Sciences  (ITS)  and  the  Australian  transmitter 
taking  part,  Backscatter  ionograms  from  the  SRI  transmitter  at  Lost  Hills,  its 
radiation  scanning  the  heated  volume  over  Platteville,  were  received/recorded  at 
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Los  Banos  while  the  Platteville  heater  operated  in  an  on-off  power  cycle.  ^ 
frequencies  ranging  from  10  to  26  MHz,  with  transmission  delays  of  about 
10  msec  corresponding  to  a  range  of  1500  km  from  WARF  sites  to  Platteville, 
the  echo  from  the  heater  induced  irregularity  volume  (heater  cloud)  is  readily 
apparent  and  quite  strong.  Figure  38  illustrates  a  nighttime  tonogram  that  shows 
a  heavy  sloping  band  of  signal  returns  from  the  ground  after  two  reflections  from 
the  ionosphere  representing  ground  backscatter,  observable  from  about  6  to 
17  MHz  over  a  two-way  transmission  delay  from  about  5  to  20  msec.  At  a  delay 
of  10  ±  0.  8  msec  the  heater  "cloud"  is  seen  as  a  "direct"  echo  at  a  target  range  of 
1500  km  for  frequencies  from  below  14  to  above  21  MHz.  For  frequencies  below 
17  MHz  the  echo  strength  of  the  heater  "cloud"  may  be  partially  suppressed  due 
to  capture  of  the  receiver's  automatic  gain  control  by  the  relatively  strong  ground 
backscatter  signal.  In  addition  to  the  "direct"  echo,  the  heater  also  produced  a 
"cloud"  echo  at  13  msec  delay  time  and  at  a  frequency  of  9  MHz  which  can  barely 
be  seen  within  the  ground  scatter.  Based  on  earlier  studies,  14  this  echo  could 
be  associated  with  a  1.  5  hop  propagation  mode.  Both  echoes  are  composed  of  two 
returns  each,  analogous  to  the  high-  and  low -angle  rays  of  oblique  propagation. 
Their  relative  strengths  depend  on  the  extent  of  defocussing  encountered. 

For  comparison,  no  effect  of  the  cloud  is  seen  during  the  heater-off  period 
at  0730  UT  (Figure  38).  The  quasi-vertical-incidence  traces  on  the  lower  left 
represent  the  ordinary-mode  trace  with  a  penetration  frequency  of  about  6.  7  MHz 
and  the  extra-ordinary  mode  trace  of  7.  3  MHz. 

3. 2. 2.  2  Australia  to  Platteville  Transmissions  (1979) 

Frequency  Modulation/Continuous  Wave  transmissions  originating  from 
Australia  were  received  at  Los  Alamos  while  the  Platteville  heater  operated  in 
an  on-off  power  cycle.  It  was  anticipated  that  the  successful  reception  of  the 
Australian  transmission  at  Los  Alamos  would  make  its  presence  over  Platteville 
likely.  There,  and  during  the  "heater-on"  part  of  the  cycle,  electromagnetic 
energy  would  be  scattered  to  the  receiving  system  at  Los  Alamos,  permitting 
deductions  about  the  prevalent  propagation  mode. 

Figure  39  illustrates  the  global  layout  of  transmitting  and  receiving  stations 
including  a  round-the-world  great -circle  path  through  Adelaide  and  Platteville. 
Portions  of  the  path  are  marked  in  2000-km  increments.  The  path  crosses  the 

13.  Basler,  RoyT.,  Showen,  Robert  L.  (1979)  Ducted  Propagation  in  the 

Ionosphere  for  ICBM  Surveillance  Application,  (U),  SRI  Int.  Tech.  Rep. 

48,  SRI  Project,  4062;  Contr.  N000  14-75-C-0930.  Funded  by 
RADC/EEP  through  ONR  NR086-876. 

14.  Bubenik,  D.  M.  (1976)  The  Combined  Effects  of  Refraction  and  Coherent 

Scattering  by  Columnar  Ionization  Density  Irregularities  in  Ionospheric 
Radio  Propagation,  Stanford  University,  Stanford,  California,  Ph.  D. 
dissertation. 
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Figure  38.  Nighttime  Ionograms  from  SRI's  Lost  Hills  Transmitter  for  27  April 
1979  (Shows  Line-of-sight  Echo  from  the  Heated  Volume  above  Platteville  When 
the  Heater  is  On  (0740  UT)  and  Off  (0735  UT) 
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Figure  39.  Map  of  the  World  Showing  the  Great- circle -path  Defined  by 
Adelaide  and  Platteville 


equator  at  a  distance  of  6800  km  from  the  location  of  the  Australian  transmitter. 
Conventional  multi -hop-propagation  between  earth  and  ionosphere  from  \ustralla 
to  Colorado  is  mostly  over  water.  Near  the  equator,  propagating  signals  are 
expected  to  encounter  the  equatorial  anomaly.  A  particular  local  mean  time  is 
indicated  on  Figure  39  corresponding  to  0700  UT.  The  day-night  terminator  is 
near  the  east  coast  of  Australia.  Most  of  the  path  to  Platteville  is  in  darkness. 

The  oblique  ionogram  in  Figure  40  illustrates  the  hop  conditions  along  the 
path  from  Australia  to  California,  including  the  occurrence  of  propagation  by  way 
of  reflection/scattering  from  the  irregularity  volume.  About  ten  modes  are  shown 
for  the  time  indicated,  with  a  spread  of  about  six  milliseconds.  When  the  heater 
Is  on,  an  additional  propagation  path  comes  into  existence  because  of  scattering 
off  the  "heater  cloud."  The  arrival  time  of  the  scattered  signal  is  delayed  by 
10  msec  relative  to  the  shortest  "direct"  path.  The  ionogram  feature  associated 
with  the  propagation  from  Australia  to  California,  by  way  of  the  irregularity 
volume  over  Platteville,  has  a  striking  similarity  to  the  corresponding  feature  of 
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Figure  40.  Oblique  Ionograms  for  the  Australia  to  Los  Banos 
Propagation  Path  With  Platteville  Heater  Off  (top)  and  On 
(bottom) 


the  heater  cloud  with  line-of-sight  illumination  observed  from  the  WARF  site 
(Figure  38).  It  was  noted  that  the  maximum  frequency  of  the  "cloud"  echo  was 
never  seen  above  the  maximum  observed  frequency  for  "direct"  propagation  at 
Los  Banos.  Moreover,  the  geometry  is  not  favorable  for  signals  coming  from 
Australia  and  arriving  horizontally  (glancing  incident)  over  Platteville,  being 
scattered  toward  California.  This  Tact,  plus  the  observed  limit  in  frequency  of 
the  echo  from  the  heater  cloud  as  seen  at  Los  Banos,  indicates  that  the  Australian 
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signal,  scattered  from  the  irregularity  volume  over  Platteville  to  Los  Banos,  did 

not  propagate  from  Australia  to  Platteville  in  an  elevated  ducted  mode,  but  rather 

13 

in  a  conventional  hop  mode.  From  Figures  35  and  36  it  may  be  surmised  that 
Los  Banos  is  indeed  not  likely  to  be  found  on  loci  of  scattering  cones  of  signals 
arriving  over  Platteville  in  elevated  ducts. 

The  lower  portion  of  Figure  41  displays  an  oblique  ionogram  for  the  trans¬ 
mission  originating  from  Adelaide,  Australia  as  received  at  Los  Alamos  over  a 
surface  distance  of  about  14,  255  km.  This  ionogram  was  obtained  on  27  April 
1979  during  the  time  interval  0800  to  0802  UT  corresponding  to  0100  to  0102  LT 
at  Los  Alamos  and  1700  to  1702  LT  at  Adelaide.  The  propagation  path  is  essen¬ 
tially  in  darkness  and  the  "terminator"  is  near  and  east  of  the  transmitter.  Sig¬ 
nals  originating  from  Adelaide  and  propagating  northeast  into  nighttime  ionosphere 
toward  Los  Alamos  and  Platteville,  encounter  negative  horizontal  ionization 
gradients  (Figure  42).  In  terms  of  ray  refraction,  such  gradients  are  favorable 
for  injection  of  a  family  or  rays  into  existing  elevated  ducts  in  which  signals 
propagate  over  large  distances  for  as  long  as  the  duct  is  effectively  maintained  in 
range,  and  with  attenuation  far  below  that  expected  for  multi-hop  propagation 
between  earth  and  ionosphere  involving  multiple  traverses  of  the  absorbing  D 
region. 

At  least  11  multi-hop  modes  are  shown  on  the  lower  left  of  Figure  41  spread 
over  a  5  msec  delay.  Relative  delay  is  read  along  the  ordinate  from  about  48  to 
55  msec.  The  time  delay  between  Salisbury  and  Los  Alamos,  resulting  from 
pulse  signal  propagating  along  the  surface  of  the  earth  over  a  great  circle  path 
with  the  velocity  of  light,  c,  in  free  space  is  obtained  dividing  the  great -circle 
distance  by  c.  The  shorter  great -circle  distance  between  Adelaide  and  Los 
Alamos  was  estimated  to  be  14,  255  km  corresponding  to  a  time  delay  of  47.  51 
msec.  With  perfect  time  synchronization  between  the  standard  clocks  at  the 
Adelaide  transmitter  and  Los  Alamos  receiver,  one  would  be  able  to  place  cor¬ 
rectly  absolute -delay  time  markers  along  the  ordinate.  This  would  allow  one  to 
determine  the  delay  in  the  arrival  of  the  earliest  earth-ionosphere  hop  mode. 

For  relatively  uniform  and  undisturbed  ionospheric  conditions  the  knowledge  of 
absolute  time  could  facilitate  mode  identification. 

The  echo  trace  extending  in  frequency  from  about  18  to  28  MHz  and  centered 
on  a  relative  delay  of  about  52  msec  throughout  this  frequency  range  with  delay 
spread  of  one  msec,  is  the  result  of  scattering  to  Los  Alamos  of  the  "Adelaide" 
signal  from  the  ionospheric  irregularity  volume  over  Platteville.  This  echo 
trace  represents  a  feature  of  those  oblique  ionograms  that  were  recorded  while 
the  Platteville  heater  was  on.  In  contrast,  the  upper  portion  of  Figure  41  illus¬ 
trates  an  ionogram  taken  5  min  earlier  when  the  Platteville  heater  was  off. 
Unrelated  to  the  heating  process,  a  change  in  the  hop-mode  structure  with  respect 
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Figure  41.  Oblique  Ionograms  for  the  Australia  to  Los  Alamos  Propagation 
Path  with  Plattevllle  Heater  Off  (top)  and  On  (bottom) 


Figure  42.  Ionospheric  Contours  Characterized  by  Plasma  Frequency  in 
Megahertz,  Derived  From  Ionospheric  Model  IONCAP  for  April,  0700  UT 
for  the  Path  from  Adelaide,  Australia  to  Los  Alamos,  New  Mexico.  A 
ray  launched  with  2  elevation  angle  at  a  frequency  of  22  MHz  is  ducted 
after  one  hop  from  the  transmitter 


to  that  shown  the  lower  portion  of  Figure  41  is  the  result  of  naturally  occurring 
ionospheric  changes  over  this  super-long  propagation  path.  The  echo  trace  pre¬ 
viously  identified  and  caused  by  heating  of  the  ionosphere  over  Platteville  is  now 
absent. 


4.  ANALYSIS  AND  INTERPRETATION 

Several  propagation  measurement  campaigns  have  been  conducted  in  con¬ 
junction  with  the  use  of  the  Platteville  heater  for  modifying  the  ionosphere,  in 
order  to  verify  the  scattering  geometry,  determine  scattering  cross  section  and 
attempt  to  extract  ducted  modes  having  propagated  over  great  distances.  These 
campaigns  provided  a  basis  for  inferring  that  ducted  modes,  scattered  by  the 
irregularity  volume,  have  been  detected  on  the  ground.  It  was  pointed  out13  that 
a  mode  scattered  from  the  irregularity  volume  over  Platteville  to  a  receiver  (at 
Los  Alamos),  where  also  regular  (direct)  multi-hop  earth-ionosphere  modes  are 
recovered,  is  likely  to  have  been  ducted  if  its  trace  on  the  obltque  ionogram  ex¬ 
tends  beyond  the  maximum  observable  frequency  of  the  multi-hop  mode  of  the 
lowest  order  (Figure  41).  If  the  maximum  observable  frequency  of  the  lowest - 
order  hop  mode  is  not  exceeded  by  the  signal  scattered  from  the  heated  volume, 
it  can  be  assumed  that  the  scattered  signal  (Figure  40)  had  arrived  at  the  heated 
volume  by  regular  hop  modes. 
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Raytracing  computations  corresponding  to  the  Australia-Platteville / 

Los  Alamos  path  were  made  at  the  Institute  for  Telecommunications  Sciences. 
Figure  42  represents  ionospheric  contours  for  April,  0700  UT,  characterized 
by  plasma  frequency  in  megahertz.  These  contours  were  provided  by  the  ITS 
program  IONCAP.  The  Australian  transmitter  is  located  at  the  origin  (D  =  0). 
Rays  are  launched  at  an  azimuth  of  65°  East  of  North  along  a  great -circle  path 
toward  Platteville.  Transmission  from  Adelaide  toward  the  East  at  0700  UT 
encounters  a  nighttime  ionosphere  with  gradients  favorable  for  the  injection  of 
rays  into  the  existing  ducts.  The  F-region  anomaly  is  seen  having  plasma  fre¬ 
quencies  of  16  MHz  at  D  =  4500  km  and  14  MHz  at  D  =  8400  km.  A  trapped  ray 
of  f  =  22  MHz  and  launch  elevation  angle  of  2°  is  illustrated.  Other  rays  were 
computed  for  26  MHz  with  14°  and  16°  elevation  angles  at  launch.  Ducted  modes 
were  obtained  in  all  three  cases.  Rays  arriving  over  Platteville  encounter 
heater -induced  field-aligned  irregularities,  that  scatter  radio  energy  to  the 
ground.  A  comparison  of  the  experimental  results  of  Figure  41  and  the  ray 
computations  of  Figure  42  provides  reasonable  certainty  that  a  ducted  propagation 
mode  was  scattered  toward  the  receiving  site  at  Los  Alamos. 

5.  CONCLUSION  , 

In  this  report  theoretical  considerations  were  presented  and  illustrated  to 
describe  expected  characteristics  of  trapped  modes  in  the  ionosphere.  Computa¬ 
tional  raytracing  through  realistic  ionosphere  models  confirmed  that  under  favor¬ 
able  gradient  conditions  and  for  certain  frequencies  and  launch  angles,  rays  can 
enter  elevated  ionospheric  ducts  and  remain  trapped  over  large  distances.  Ex¬ 
perimental  evidence  of  scattering  from  artificially  created  ionospheric  irregular¬ 
ities  showed  that  modes  trapped  in  such  elevated  ionospheric  channels  can  be 
detected  on  the  ground  at  locations  that  generally  satisfy  the  geometry  for  scat¬ 
tering  from  field-aligned  irregularities. 

While  these  ground-based  measurements  have  provided  new  understanding 
and  evidence  for  ionospheric  ducting,  the  experimental  configuration  was  fixed 
and  thus  limited.  A  new  approach  is  presently  envisioned  and  plans  have  been 
Implemented  to  conduct  a  ducting  experiment  using  a  ground-based  transmitter 
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with  the  receiver,  mounted  on  a  space  vehicle,  orbiting  at  altitudes  of  the  ducting 
channel.  It  will  then  be  possible,  for  the  first  time,  to  measure  in  situ  the  extent 
of  ducting  in  ionospheric  channels,  and  the  properties  of  ionospherically  ducted 
modes. 
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